Many-body effects in transport through open systems: pinning of resonant levels 
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The role of electron-electron interaction in transport properties of open quantum dots is stud- 
ied. The self-consistent full quantum mechanical magnetotransport calculations within the Hartree, 
Density Functional Theory and Thomas-Fermi approximations were performed where a whole de- 
vice, including the semi-infinitive leads, is treated on the same footing (i.e. the electron-electron 
interaction is accounted for both in the leads as well as in the dot region). The main finding of the 
present paper is the effect of pinning of the resonant levels to the Fermi energy due to the enhanced 
screening. Our results represent a significant departure from a conventional picture where a varia- 
tion of external parameters (such as a gate voltage, magnetic field, etc.) causes the successive dot 
states to sweep past the Fermi level in a linear fashion. We instead demonstrate highly nonlinear 
behavior of the resonant levels in the vicinity of the Fermi energy. The pinning of the resonant levels 
in open quantum dots leads to the broadening of the conduction oscillations in comparison to the 
one electron picture. The effect of pinning becomes much more pronounced in the presence of the 
perpendicular magnetic field. This can be attributed to the enhanced screening efficiency because of 
the increased localization of the wave function. The strong pinning of the resonant energy levels in 
the presence of magnetic field can have a profound effect on transport properties of various devices 
operating in the edge state transport regime. We also critically examine an approximation often 
used in transport calculations where an inherently open system is replaced by a corresponding closed 
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I. INTRODUCTION 

A transport regime where a sub-micron lateral struc- 
ture is strongly coupled to electron reservoirs (leads) is 
usually referred to as an open onei^. This transport 
regime can be realized in quantum wires, dot and anti- 
dot structures typically fabricated using a split-gate or 
related techniques. Such the techniques allow one to ob- 
tain devices with desired and variable geometry and pa- 
rameters such as an electron density and lead openings. 
The quantum dot operates in an open regime when the 
gate voltage sets up two quantum point contacts (QPCs) 
at the entrance and exit of the dot such that they trans- 
mit one or more channels (i.e. the conductance of an 

individual QPC, Gqpc > tt)- I n this regime electrons 
can freely enter and exit the dot, such that the electron 
number inside the dot is not integer and the chemical po- 
tential throughout the whole device in the linear responcc 
regime is constant. [Opposite regime emerges when the 
point contacts are nearly pinched off and is referred to as 
a Coulomb blockade^. In this case the electron number 
in the dot is quantized and the chemical potential inside 
the dot is different from that one in the leads]. During 
the last decade the open quantum dots have received a 
significant attention providing many important insights 
into areas such as quantum interference, chaos, decoher- 
ence, localization and many others^. Earlier transport 



experiments have been mainly devoted to the dots with 
hundreds or even thousands electrons. Only recently it 
has become possible to reduce occupancy down to only 
a few or even one electronic. 

Electron-electron interaction is known to have great 
impact on transport in quantum dots with such pro- 
nounced examples as Coulomb blockaded or Kondo 
effect^. A description of the quantum transport in quan- 
tum dots is often based on model Hamiltonians con- 
taining phcnomenological parameters such as coupling 
strengths or charging constants^^. In many cases it 
is not always straightforward to relate quantitatively the 
above parameters to the physical processes they represent 
in the real system and sometimes it is not even obvious 
whether a model description is sufficient to capture the 
essential physics. At the same time, it is now becomes 
well recognized that a detailed understanding and inter- 
pretation of the experiment might require a quantitative 
microscopical modelling of the system at hand free from 
phenomcnological parameters and not relying on model 
Hamiltonians which validity is poorly controlled. The 
importance of such the modelling can be illustrated by 
examples including the quantitative description of the 
compressible/incompressible strips in magnetic field at 
the edges of the two-dimensional electron gas^ or expla- 
nation of the Hund rule observed in fcw-clcctron quantum 
dotsii, just to name a few. 
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The purpose of the present paper is twofold. First, we 
develop an approach aimed on full quantum mechanical 
many-body transport calculations in open systems that 
starts from the lithographical layout of the device and 
does not include phcnomenological parameters such as 
coupling strengths, charging constants etc. The whole 
device, including semi-infinitive leads, is treated on the 
same footing (i.e. the electron-electron interaction is ac- 
counted for both in the leads as well as in the dot region) . 
Using the recursive Green's function technique we self- 
consistently compute the scattering solutions of the two- 
dimensional Schrodinger equation in magnetic fielcU 2 -. 
Following the parametrization for the exchange and cor- 
relation energy functionals of Tanatar and Ceperley^ 
the electron-electron interaction is incorporated within 
the density functional theory (DFT) in the local density 
approximatio n 11 ' 1 . The validity of the DFT approxi- 
mation is supported by the excellent agreement with the 
exact diagonalization and variational Monte Carlo cal- 
culations performed for few-electron systems^, as well 
as by the good quantitative correspondence between the 
experiment and the DFT calculations for the magneto- 
conductance of quantum wires^. Our approach thus ac- 
counts for the quantum mechanics nature of the scat- 
tering states and the resonant levels as well as for the 
exchange and correlation beyond the Hartree approxi- 
mation. At the same time, it accurately describes the 
global electrostatic and the screening in the dots as well 
as in the leads. 

The second aim of the present paper is to revise the 
role of the electron-electron interaction in transport prop- 
erties of open quantum dots. It is widely believed that 
in open transport regime (as opposed to the Coulomb 
blockade or Kondo regime), the electron-electron inter- 
action plays only a minor role. The main finding of the 
present paper is the effect of pinning of the resonant lev- 
els to the Fermi energy due to the enhanced screening. 
Our results represent a significant departure from a con- 
ventional picture adopted in most model Hamiltonians 
as well as in more sophisticated numerical calculations 
where a variation of external parameters (such as a gate 
voltage, magnetic field, etc.) causes the successive dot 
states to sweep past the Fermi level in a linear fashion. 
We instead demonstrate highly nonlinear behavior of the 
resonant levels in the vicinity of the Fermi energy. One 
of the observable consequence of this effect is smearing 
of the conductance fluctuations. We also show that the 
resonant level pinning becomes especially pronounced in 
magnetic field. Thus, accounting for this effect might 
be important for the interpretation of the magnctotrans- 
port experiment in open structures, including e.g. recent 
studies of the electronic Mach-Zehnder interferometer— 
and the Laughlin quasiparticle interferometer—, struc- 
tures designed to test the realization of the topological 
quantum computing^, antidot struct ure a 20 ' 21 ' 22 and oth- 
ers. It should be also noted that the quantum dot struc- 
tures demonstrating Kondo effect fall into the semi-open 
transport regime such that accounting for effect of the 



nonlinear screening leading to the resonant level pinning 
might be essential for the interpretation of the experi- 
ments in this regime as well. 

The paper is organized as follow. Section II presents 
the model and the Hamiltonian of the system at hand. In 
Sec. Ill the numerical method for the self-consistent cal- 
culation of the magnctoconductance is described. Com- 
putational results are presented and discussed in Sec. IV, 
and conclusions are given in Sec. V. 

II. MODEL 

We consider an open quantum dot attached to semi- 
infinitive leads (electron reservoirs) in a perpendicular 
magnetic field B. A schematic layout of the device is il- 
lustrated in Fig. [lja). Charge carriers originating from a 
fully ionized donor layer form the two-dimensional elec- 
tron gas (2DEG) which is buried inside a substrate at the 
GaAs/Al^Gai-^As hctcrointcrface situated at the dis- 
tance b from the surface. Metallic gates placed on the 
top of the heterostructure define the dot and the leads 
on the depth of the 2DEG (Figs. [1] (a),(b)). 

The Hamiltonian of the whole system (the dot + the 
leads) can be written in the form 

H = H + V(r) (1) 

where Hq is the kinetic energy in the Landau gauge, A = 
(-By, 0,0), 

fa? | f d ciBy \ ^ d^ 1 1 

r = (x,y), m* = 0.067TO e is the GaAs effective mass. 
The total confining potential within the framework of 
the density functional theory is the sum of the electro- 
static confinement potential, the Hartree potential, and 
the exchange-correlation potential 

V(r)=V conf (r)+V H (r)+V xc (r)- (3) 

The electrostatic confinement V co „/(r) = Vg a tes{^) + 
V d onors + Vschottky includes contributions respectively 
from the top gates, the donor layer and the Schottky bar- 
rier. The explicit expressions for the potentials V 9ates (r) 
and Vdonors are given in respectively Ref. and Ref. 
l24l ; the Schottky barrier is chosen to be Vschottky = 0.8 
cV. The Hartree potential is written in a standard form 

Vff(r) = — — [ dr'n(r') ( - 1 = ) , 

(4) 

where n(r) is the electron density and the second term 
describes the mirror charges places at the distance of b 
from the surface, n(r) is the electron density, e r = 12.9 
is the dielectric constant of GaAs, and the integration is 
performed over the whole device area including the semi- 
infinite leads. 
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FIG. 1: (color online) (a) Structure of an open quantum dot. 
The internal region is attached to two semi-infinitive quantum 
wires which serve as the electron reservoirs, (b) The layout of 
the gates defining the dot. The dashed line defines the area S 
of the quantum dot used to calculate a number of electrons in 
it. (c) Representative calculated the self-consistent potential 
and the electron density. 



The last term in the total confining potential Q is the 
exchange and correlation potential T4 c [?i(r)] = V^[n(r)] + 
V c [n(r)] which is the functional of the electron density. In 
the local density approximation it is given by a functional 
derivative^ 



V xc = -j- {nt xc (n)} . 



(5) 



For e xc we have used the parametrization of Tanatar and 
Ceperleyi^ for the case of spin degenerate electrons. In 
particular, for the exchange potential this parametriza- 
tion gives 14[n(r)] 
ting V xc (r) 



- eo£ e 7r3/2 y / n(r)/2. Note that sct- 
in Eq. ([3]) we reduce our approach to the 



standard Hartrec approximation. 

To outline the role of quantum mechanical effects in 
the electron-electron interaction in open quantum dots 



we also consider the Thomas-Fermi (TF) approximation. 
In this approximation the kinetic energy is related to the 
electron density^, H n = ^-n(r). The self-consistent 
electron density is thus obtained from the solution of the 
equation 



nh 2 

—n(x,y) + V conf (r) + V H (r) 



E F . 



(6) 



The electron density and the total confining potential 
calculated within the TF approximation do not capture 
quantum- mechanical quantization of the electron motion. 
The utilization of the TF approximation for the mod- 
elling of the magnctotransport in open system is therefore 
conceptually equivalent to a one-electron approach. The 
difference between these approaches is the shape of the 
total confining potential: in one-electron transport sim- 
ulations one typically starts with a model hard-wall con- 
finement, whereas the TF approximation gives a rather 
smooth potential which represents a good approximation 
to the actual confinement. 



III. METHOD 

The magnetoconductance through the quantum dot 
in the linear response regime is given by the Landauer 
formula— 



E F ) 



(7) 



where T(E) is the total transmission coefficient, fp d(E — 
Ep) is the Fermi-Dirac distribution function and Ep 
is the Fermy energy. In order to calculate T(E) in 
perpendicular magnetic field we utilize the recursive 
Greens function technique in the hybrid energy-space 
representation— We discretize Eq. §T§ and introduce the 
tight-binding Hamiltonian (with lattice constant a = 4 
nm), where the perpendicular magnetic field is included 
in a form of the Pcicrl's substitution—. The retarded 
Green's function is introduced in a standard way2£, 



G = (E-H + irjY 



(8) 



The Green's function in the real space representation, 
Q(r,r, E), provides an information about the electron 
density at the site 



n(r) = --3 J dEg(r,r,E) f(E - Ep). 



(9) 



Note that Q(r,r,E) is a rapidly varying function of en- 
ergy. As a result, a direct integration along the real axis 
in Eq. ^ is rather ineffective as its numerical accuracy is 
not sufficient to achieve convergence of the self-consistent 
electron density. Because of this, we transform the inte- 
gration contour into the complex plane 5s[E] > 0, where 
the Green's function is much more smoother. Note that 
all poles of the retarded Green's function are in the lower 
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FIG. 2: A typical integration contour used in the calculation 
of integral Dots indicate the poles of the Fermi-Dirac 

distribution function in the upper complex plane at J? [E] = 
E F , 9 [E] = (2m + l)TvkT, m = 0, 1, 2, .... V b is the bottom 
of the conduction band in the leads (the lowest potential in 
the system). 



half-plane ^s[E] < 0. A typical contour used in the in- 
tegration avoiding poles of the Fermi-Dirac function is 
shown in Fig. [2] 

In order to calculate the Green's function of the whole 
system (dot + leads) we divide the system into three 
parts, the internal computational region and two semi- 
infinite leads as shown in Fig. [1] (a). Note that the in- 
ternal region incorporates not only the dot, but also the 
straight segments including a part of the leads. We place 
the semi-infinite leads sufficiently far away from the dot 
where the total self-consistent potential and the electron 
density do not change along the leads (i.e. the electron 
density and the potential in the leads are not affected by 
the internal region, such that the leads can be considered 
as uniform quantum wires.). This allows us to calculate 
the total potential and the electron density in the lead 
regions using the technique developed in Ref. [26| for an 
infinite homogeneous channel in the perpendicular mag- 
netic field. The Green's function in the internal region is 
calculated using the standard recursive Green's function 
technique. The total Green's function for the whole sys- 
tem is calculated by linking with the help of the Dyson 
equation the surface Greens function for the semi-infinite 
leads (with the self-consistent potential calculated using 
the technique of Ref. [2fjh and the Green's function of the 
internal region. 

All the calculations described above are performed self- 
consistently in an iterative way until a converged solution 
for the electron density and potential (and hence for the 
total Green's function) is obtained. Having calculated 
the total self-consistent Green's function, the scattering 
problem is solved where the scattering states in the leads 
(both propagating and evanescent) are obtained using 
the Green's function technique of Ref. [26| (The equation 
for the calculation of the transmission and reflection co- 
efficients using the Green's function in the presence of 
the magnetic field is derived in Ref. [r3 ). 

Having calculated the Green's function of the internal 
region and the wave function in the leads we can recover 
the wave function ?/>(r, E) inside the internal region. For 
visualization of the wave function inside the dot we in- 



clude the effect of the finite temperature as follows, 

\Hr)\ 2 = -JdE\^Ef dfFD ^ E EF) . (10) 

In order to find the DOS inside the quantum doli 2 ^ we 
perform integration over the dot area S defined in Fig. 
HTb), 



DOS(E) = --$ / drg(r,r,E). 
t Js 



(11) 



Note that a many-body approach conceptually similar 
to ours for calculation of the quantum transport in an 
open dot was developed in Ref. H3 for the case of zero 
magnetic field. The magnetic field was included in the 
dot region in the transport calculations presented in Ref. 
l28l where, however, the leads were considered as non- 
interacting. 

As we mentioned in the previous section, we also em- 
ploy the TF approximation for the calculation of the 
magnetotransport through the quantum dot. This cal- 
culation is done by the same method as described above, 
with the only difference that the self-consistent electron 
densities in the internal region and in the semi-infinite 
leads are calculated from the semi-classical FT equation 
©, as opposed to Eq. ^ that relates the electron den- 
sity to the quantum mechanical Green's function. 

The self-consistent solution in quantum transport or 
electronic structure calculations is often found using a 
"simple mixing" method, when on each m + 1 iteration 
step a small part of a new potential is mixed with the 
old one (from the previous iteration step), V m +i = (1 — 
e)V m + e V n+1 , e being a small constant, ~ 0.1 — 0.01. It 
is typically needed ~ 200 — 2000 iteration steps to achieve 
our convergence criterium 



m+l 



,m+ 1 



< 10" 



(12) 



where n m = J n™ (r) dy is the electron density in the 
dot on m-th iteration step. In order to improve the 
convergence we employ the modified Broyden's second 
method^ which allows us to reduce drastically the num- 
ber of iteration steps to ~ 15 — 40. An input charge den- 
sity for m + 1 iteration is constructed from the 
sets of input and output densities (n™ t , n™), from all m 
previous iterations 



nZ +1 



B ipm _ u j (V j ) F m (13) 
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The initial guess B l is taken to be a small constant 
(~ 0.1 — 0.01) so that the input to the second iteration 
is effectively constructed using the simple mixing. The 
scheme given by Eqs. (fT^j) requires the storage of rela- 
tively small number vectors which is thus more effective 
than the original Broyden's second method. 

IV. RESULTS AND DISCUSSION 
A. A few electron open dot 

We calculate the magnetotransport of a split-gate open 
quantum with following parameters representative for a 
typical experimental structure. The 2DEG is buried at 
b = 60 nm below the surface (the widths of the cap, 
donor and spacer layers are 10 nm, 36 nm and 14 nm 
respectively), the donor concentration is 0.6 • 10 24 m~ 3 . 
The width of the semi-infinitive leads is wi eac i = 540 nm, 
and the width of the constrictions is w qpc = 100 nm (both 
the quantum point contacts are identical), see Fig. [1] 
(b). The length of the quantum dot is kept constant 
throughout the paper, Idot = 160 nm, while the width of 
the dot is varied in the range Wdot = 170 — 440 nm. The 
gate voltages applied to the gates is Vi ea d = —0.4 V and 
V q pc = —0.44 V. With these parameters of the device 
there are 14 channels available for propagation in the 
leads and the electron density in the center of the leads 
is ni ea d = 1.6 • 10 15 m~ 2 . The maximal electron density 
in the dot (for Wdot = 440 nm) is also ridot = 1-6 • 10 15 
m~ 2 . The temperature is fixed at T = 0.2 K for all results 
presented below. 

In the following we discuss the open quantum dot with 
N = 1 propagating channel through both quantum point 
contacts. Note, that increasing the number of propagat- 
ing channels to TV = 2,3 does not qualitatively change 
the results presented below. In order to set up the QPCs 
in the one-mode regime, we grounded one of them and 
studied the conductance of the remaining QPC as a func- 
tion of the gate voltage V qpc . The calculated conduc- 
tance shows a characteristic step-like dependence and 
we choose V qpc at the first conductance plateau, namely, 
V qpc = -0.44 V. 

Figure [3] shows the color-scale plot of the conductance 
G as a function of the dot width Wdot and the gate volt- 
age V g (the magnetic field is restricted to zero). The 
distinctive feature of the open quantum dots is the os- 
cillations of the conductance, in response to the change 
of the geometrical size or the Fermi energy. (Various as- 
pects of the conductance oscillations in small and large 
dots have been the subject of numerous experimental and 
theoretical works during the last decade 1 -) . In the present 
study we concentrate on two dots with Wdot = 170 nm 
and Wdot = 430 nm (as indicated by arrows in Fig. [3]). 
As will be shown below, the first dot corresponds to the 
transport regime when A 3> V, whereas the second dot 
operates in the regime A ~ T, where A = ^Vg is the 
mean level spacing separation in the dot with the actual 
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FIG. 3: (color online) The conductance of the open quantum 
dot as a function of the width, Wdot, and the gate voltage V g 
calculated within the Hartree approximation. Solid thin lines 
denote the number of electrons. [Note that a zig-zag-type 
behavior of the electron number is an artifact due to finite grid 
steps]. Arrows indicate Wdot = 170 nm and Wdot = 430 nm 
corresponding to two regimes, A>F and A ~ V, discussed 
in Sec. IV. 



area S a , and T is the lead- induced broadening of the res- 
onant energy levels. 

It is worth to note that a linear change of the gate volt- 
age leads to a nonlinear change of the effective dot size 
(If it were linear, the conductance oscillations in Fig. [3] 
would exhibit a perfect linear-stripe-type pattern). This 
is consistent with the experimental results of Ref. [30 that 
show a deviation from the linear dependence of the gate 
depletion distance as the gate voltage was varied. 



B. Regime A»T 

Figure [4] shows the number of electrons in the few- 
electron open quantum dot, the conductance and the 
peak-energy level position as a function of the gate volt- 
age V g for a quantum dot with Wdot = 170 nm calcu- 
lated in the Hartree, DFT and TF approximations (The 
Hartree approximation corresponds to disregarding of the 
exchange-correlation potential in Eq. V xc = 0). The 
peak positions of the resonant energy levels in the dot are 
extracted from the calculated DOS, as illustrated in Fig. 
131(c) for the case when the gate voltage V g = —0.6 V. The 
estimation of the mean level separation gives A w 0.6 
meV (for V g = —0.6 V) which agrees quite well with the 
actual level separation shown in Fig. [4] (c) - (e). [Note 
that within the given interval of variation of V g , the ac- 
tual dot area changes and hence A varies as well]. An 
inspection of the DOS shows that the the separation be- 
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FIG. 4: (color online) (a) The number of electrons, (b) the 
conductance, and (c)-(e) resonant energy structure in the 
few-electron open quantum dot with Wd o t = 170 nm calcu- 
lated within the Hartree, DFT and Thomas-Fermi approxi- 
mations. The top panel shows the electron probability am- 
plitudes |\P(a;,y)| 2 for the resonant energy levels marked by 
arrows in (c). Inset in (c) shows the DOS for V g —-0.6 V. 



tween the resonant levels are much larger than the level 
broadening, A>T. 

All three approximations give very similar electron 
number N in the dot as a function of the gate voltage 
V g . N monotonically increases with increase of V g , which 



reflects the fact that electrons freely enter and leave the 
dot in the open regime. However, the conductance calcu- 
lated in the Hartree and the DFT approximations exhibit 
rather similar behavior, whereas the Thomas-Fermi con- 
ductance shows a very different gate voltage dependence. 
The origin of this difference can be understood from the 
analysis of the resonant level structure. The resonant 
energy levels calculated within the Hartree and DFT ap- 
proximations get pinned to the Fermi energy, see Figs. 
SJc) and (d) respectively. In stark contrast, the energy 
level positions calculated in the Thomas-Fermi approx- 
imation sweep past the Fermi level in a linear fashion 
when the applied voltage is varied, Fig. [4] (e). 

The effect of level pinning is related to the screening 
properties of the open quantum dot and the presence of 
the resonant structure in its DOS. Indeed, in the vicin- 
ity of the resonances the DOS of the dot is enhanced 
such that electrons with the energies close to Ep (when 
/fd < 1) can easily screen the external potential. This 
leads to the metallic behavior of the system when the 
electron density in the dot can be easily redistributed to 
keep the potential constant. As a result, in the vicinity 
of a resonance the system only weakly responds to the 
external perturbation (change of a gate voltage, mag- 
netic field, etc.), i.e. the resonant levels becomes pinned 
to the Fermy energy. A comparison between the Hartree 
and the DFT approximations indicate that the exchange- 
correlation interaction seems to enhance the pinning, but 
an overall change is small (c.f. Figs. 0](c) and (d)). Thus, 
in the following we concentrate on the Hartree approx- 
imation only. In contrast to the Hartree and the DFT 
approaches, the effect of pinning is absent in the TF ap- 
proximation. This is because the total confining potential 
calculated within the TF approximation does not capture 
the resonant structure of the DOS. Note that a modelling 
of the magnetotransport in a quantum dot conceptually 
similar to our TF approach was performed by Bird et 
al£^, where the confining potential for every given gate 
voltage was obtained as a self-consistent solution of the 
Poisson equation. 

Conductance oscillations in the open quantum dots can 
be related to presence/absence of the resonant energy lev- 
els at Ep. Indeed, the first three peaks in the Hartree 
conductance are attributed to the presence of correspond- 
ing resonant levels, which is confirmed by the inspection 
of the wave function |\I/(a;,y)| 2 (by counting the number 
of nodes of |^(a;,j/)| ), Fig. 2J In turn, the dips indi- 
cate the absence of levels at Ep and agree precisely with 
integer N reflecting the fact that all the available lev- 
els are below Ep and thus are fully filled. Note that for 
less negative gate voltages (|V| < 0.3 V) the separation 
between the levels A becomes comparable to the broad- 
ening r, such that the dips in the conductance are no 
longer correspond to the integer electron number TV (see 
next section for details) . It is also interesting to note that 
the resonant state corresponding to the 5th eigenstate (5 
nodes of the \^(x, y)\ 2 ) is situated lower in energy than 
the corresponding resonant state related to the 4th eigen- 
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state. This is an indication that 5th state couples with 
the leads more strongly that 4th stated. 

The pinning of the resonant energy levels has an im- 
portant effect on transport in open quantum dots. In the 
considered transport regime, A 3> T, the conductance 
calculated within the one-electron (TF) approximation 
exhibits distinct peaks separated by broad valleys of es- 
sentially zero conductance. This reflects the structure 
of the one-electron DOS where the resonant levels sweep 
past the Fermi level in a linear fashion. In contrast, as 
a result of pinning, the DFT and Hartree conductance 
shows much broader oscillations in comparison the one- 
electron approach (Fig. Bib)). 

C. Regime A ~ V 

When an effective size of a quantum dot increases, the 
mean level spacing separation A decreases. For the quan- 
tum dot of the size Wdot = 430 the estimation of the mean 
level separation gives A ss 0.15 meV (for V g = —0.6 V) 
which agrees quite well with the actual level separation 
shown in Fig. (c) - (d). An inspection of the DOS 
shows that for this dot the spacing between neighboring 
levels is comparable with the level broadening, A ~ T, 
see Fig. [SJc). [Note, that the level broadening is con- 
trolled by coupling to the leads and does not depend on 
the dot size.] Because neighboring levels start to over- 
lap, there is always one or several energy states at Ep 
mediating transport through the dot. The effect of the 
pinning of the energy levels in the Hartree calculations 
(as well as its absence in the TF calculations) is also 
clearly seen in this regime as well. However, in contrast 
to the regime A>T, the dips in the conductance can no 
longer be related to the integer number of electrons in the 
dot. Instead, it is the interference between states at the 
entrance and exit of the open quantum that determines 
the dependence G = G(V g )^. Our preliminary results 
for the conductance of even larger dots (containing hun- 
dreds of electrons) outline different features of the TF 
and Hartree conductances that are manifest themselves 
in the amplitude and the broadening of the conductance 
peaks. A detailed analysis of the statistics of the con- 
ductance oscillations is however outside the scope of the 
present work and will be deferred to future publications. 




V (V) 

s s 

FIG. 5: (color online) (a) The number of electrons, (b) the 
conductance, and (c)-(d) resonant energy structure in the 
few-electron open quantum dot with Wdot =430 nm calculated 
within the Hartree and Thomas-Fermi approximations. The 
electron probability amplitudes [^(x, y)\ 2 are shown for some 
representative V g (top panel). Inset in (c) shows the DOS for 
V 9 =-0.5 V. 



D. Effect of magnetic field 

In a sufficiently high magnetic field the electron trans- 
port takes place by the edge states with a characteris- 
tic dimension of the order of the magnetic length Ib = 
^Jh/eB. In the edge state transport regime backscatter- 
ing on the potential defining the quantum dot decreases 
and, for a large enough B 7 electrons pass through the 
device with the transmission close to unity. Transport 
in such a regime is referred to as adiabatic. For the 
open quantum dot of Wdot = 430 nm, transition to adi- 



abatic propagation takes place at about B » 0.5 T, see 
Fig. [5{a). The conductance for B > 0.5 T shows pro- 
nounced oscillations due to the Aharonov-Bohm inter- 
ference. When the magnetic field changes such that the 
total magnetic flux $ = BS through the dot modifies by 
the one flux quantum (j)Q = h/e, the conductance demon- 
strates periodic oscillations with the period AB = <f>o/S 
(S is the characteristic area of the dot). Using the actual 
dot area S a we get AB = 0.11 T, which is nearly twice 
less than extracted from Fig. E^a), where AB = 0.25 
T. The discrepancy can be related to a finite extent of 
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FIG. 6: (color online) (a) The conductance of the open quan- 
tum dot as a function of the magnetic field B and gate voltage 
V g for the open quantum dot of the width u>dot=430 nm cal- 
culated within the Hartree approximation, (b) The energy 
structure for B = 0.5 T calculated within the Hartree ap- 
proximation. The top panel shows the electron probability 
amplitudes \^{x,y)\ 2 (top) for two representative magnetic 
fields. 



the edge state circulating inside the dot (Ib ~ 35 nm for 
B = 0.5). As a result, the area enclosed by the edge state 
is much smaller than the geometrical area of the dot. 

The resonant energy structure is modified substan- 
tially when a magnetic field is applied, c.f. Fig. [6](b) 
for B = 0.5 T and Fig. [5] (c) for B = T. The res- 
onant levels exhibit almost equal separation which can 
be related to the well-known Darvin-Fock-type energy 
spectrum formation for the corresponding closed dofc22 
The distinguished feature of the energy level structure 
is much stronger pinning of the resonant levels to Ep 
that persists over larger intervals of V g in comparison 
to the B = case. Stronger pinning can be attributed 
to the enhanced screening efficiency because of the in- 
creased localization of the wave function for the case of 
nonzero magnetic field. As we mentioned in the introduc- 
tion the strong pinning of the resonant energy levels in 
the presence of the magnetic field can have a profound ef- 
fect on thransport properties of various devices operating 



in the edge state transport regime including the Mach- 
Zender— and the Laughlin^ interferometers as well as 
antidot device o 19 i 20 ' 21 i 22 . 

To conclude this section, it is worth mentioning that 
another manifestation of the screening in the edge state 
regime is the well-known effect of formation of the com- 
pressible and incompressible strips near the structure 
boundary^. 

E. Open vs closed system 

In modelling quantum mechanical transport in quan- 
tum dots and related systems one often uses an approx- 
imation where an inherently open system is replaces by 
a corresponding large, but closed one, see e.gi^i. In this 
section we critically examine such an approximation. In 
particular, we address a question whether a conductance 
calculated in such a way coincides with the conductance 
of a truly open system, and whether the pinning of the 
resonant energy levels is survived or not. For modelling 
of the closed system, we replace the semi-infinitive leads 
by the potential walls of an infinitive height such that 
the solution of the Schrodinger equation reduces to the 
eigenproblcm 

Hij) = Ei/>, (14) 

where E and tp arc discrete sets of eigenvalues and eigen- 
vectors, and the Hamiltonian H is given by Eq. (JTJ) . 
We solve Eq. (fT4|) self-consistently performing the fast 
Fourier transformation from the coordinate into the en- 
ergy space, which greatly reduces a computational cost. 

In our calculations we fix the Fermi energy such that 
the charge density in the system is given as 

n{v)=Y J ^)f{E i -E F ). (15) 

i 

This assumption leads to a noninteger electron number 
in a system, but we apriory construct the closed system 
resembling the open one as much as possible. Note also, 
that because the total number of electrons N to t 3> 1, the 
effect of the non-intcger N to t on the total potential is 
practically negligible. 

In order to calculate the conductance of the system at 
hand, we cut off slices in the vicinity of the boundaries 
as illustrated in Fig. [T^a) and then add homogeneous 
semi-infinite leads with the potential that matches the 
potential of the boundary slices. Finally, we solve a scat- 
tering problem for this given potential using the recursive 
Greens function technique^ 2 -. 

Figures E)Jb)-(d) show the electron number in the dot 
N, the conductance G and the resonant energy struc- 
ture within the Hartree approximation. The comparison 
to the corresponding results for the open system shows 
that the closed-system approximation reproduces all the 
results not only qualitatively but rather quantitatively, 
c.f. Fig. 2]Thc on ly difference is the shift along V^-axis 
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FIG. 7: (color online) (a) The representative self-consistent 
charge densities for the closed and open systems. The thick 
dashed lines show the cuts for the transport calculations, (b) 
The electron number, (c) the conductance and (d) the peak 
energy-level position calculated within the Hartree approxi- 
mation in the closed-system approximation. The dashed lines 
denote the result for the open system (the same as in Fig. [3} . 
The width of the quantum dot is Wdot~170 nm. 



which is simply related to the fact that the Hartree po- 
tential for the case of the closed system, in contrast to the 
open one, does not include a contribution from the semi- 
infinite leads. Figure[7]Jd) reveals that the pinning of res- 
onant energy levels to the Fermi energy is present in the 
closed system-approximation as well. We thus conclude 
that this approximation might be used for the modelling 
of the transport properties and the resonant energy level 
structure of the corresponding open system. We how- 
ever should note that with the present approximation we 
could not satisfy the convergence criterium (fT2"|) that we 
routinely use for the calculation of the conductance in 
the open systems as described in Sec. III. 



V. CONCLUSION 

We have developed an approach for full quantum me- 
chanical many-body magnctotransport calculations in 
open systems that starts from the lithographical layout of 
the device and does not include phenomenological param- 
eters like coupling strengths, charging constants etc. The 
whole device, including semi-infinitive leads, is treated on 
the same footing (i.e. the electron-electron interaction is 
accounted for in both the leads as well as in the dot re- 
gion). The many-body effects are included within the 
DFT and the Hartree approximations. 

The developed method was applied to calculate the 
conductance through an open quantum dot. The main 
finding of the present paper is the effect of pinning of the 
resonant levels to the Fermi energy due to the enhanced 
screening. Our results represent a significant departure 
from a conventional picture where a variation of external 
parameters (such as a gate voltage, magnetic field, etc.) 
causes the successive dot states to sweep past the Fermi 
level in a linear fashion. We instead demonstrate highly 
nonlinear behavior of the resonant levels in the vicinity 
of the Fermi energy. We show that the pinning effect is 
absent in a one-electron (Thomas-Fcrmy) approximation 
because in this case the self-consistent potential does not 
account for the resonant structure of the DOS in the dot. 
The pinning of the resonant levels in open quantum dots 
leads to the broadening of the conduction oscillations in 
comparison to the one electron picture. It remains to be 
seen whether accounting for this effect might shed new 
light on the interpretation of the conductance oscillation 
statistics in open quantum dots. 

The pinning of the resonant levels becomes much more 
pronounced in the presence of the perpendicular mag- 
netic field. This can be attributed to the enhanced 
screening efficiency because of the increased localization 
of the wave function. The strong pinning of the reso- 
nant energy levels in the presence of magnetic field can 
have a profound effect on transport properties of vari- 
ous devices operating in the edge state transport regime 
including Mach-Zender— and Laughlir>i£ interferometers 
as well as antidot device a 19 i 20 ' 21 i 22 . 

We should stress that the pinning effect predicted in 
this paper is not specific to the considered material sys- 
tem (GaAs/AlGaAs heterostructure) and is expected to 
hold in any two-dimensional system in open transport 
regime (e.g. Si inversion layer structures, etc.). 

Finally, in the present paper we critically examined an 
approximation used in modelling of quantum mechanical 
transport in quantum dots and related systems when an 
inherently open system is replaces by a corresponding 
large, but closed one. 

In the present study we have limited ourselves to the 
case of spinless electrons. Work is in progress to include 
the effect of the spin in order to revisit the effect of spin 
splitting recently observed in open quantum dots 9 -^. 
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